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The ground-state properties of superfluid nuclear systems with 1 So pairing are 
studied within a local energy-density functional (LEDF) approach. A new form of 

^\J ' the LEDF is proposed with a volume part which fits the Friedman-Pandharipande 

and Wiringa-Fiks-Fabrocini equation of state at low and moderate densities and 
allows an extrapolation to higher densities preserving causality. For inhomoge- 
neous systems, a surface term is added, with two free parameters, which has a 
^, ■ fractional form like a Pade approximant containing the gradient of density squared 

in both the numerator and denominator. In addition to the Coulomb direct and 
exchange interaction energy, an effective density-dependent Coulomb-nuclear cor- 
relation term is included with one more free parameter. A three-parameter fit 
to the masses and radii of about 100 spherical nuclei has shown that the latter 
term gives a contribution of the same order of magnitude as the Nolen-Schiffer 
anomaly in Coulomb displacement energy. The root-mean-square deviations from 
experimental masses and radii with the proposed LEDF come out about a fac- 
tor of two smaller than those obtained with the conventional functionals based on 
the Skyrme or finite-range Gogny force, or on the relativistic mean-field theory. 
The generalized variational principle is formulated leading to the self-consistent 
Gor'kov equations which are solved exactly, with physical boundary conditions 
both for the bound and scattering states. The method is used to calculate the 
differential observables such as odd-even mass differences and staggering in charge 
(-H | radii. With a zero-range density-dependent cutoff pairing interaction incorporating 

a density-gradient term, the evolution of these observables is reproduced reason- 
ably well, including kinks at magic neutron numbers and sizes of staggering. An 
extrapolation from the pairing properties of finite nuclei to pairing in infinite nu- 

^%J ' clear matter is discussed. A "reference" value of the pairing gap Ap ss 3.3 MeV 

is found for subsaturated nuclear matter at about 0.65 of the equilibrium density. 
With the formulated LEDF approach, we study also the dilute limit in both the 
weak and strong coupling regime. Within the sum rules approach it is shown that 
the density-dependent pairing may also induce sizeable staggering and kinks in the 
evolution of the mean energies of multipole excitations. 

1 Introduction 

The challenge of deriving the properties of nuclear matter and finite nuclei starting 
with a "realistic" bare NN interaction stimulates broad activity in developing the 
microscopic approaches to strongly coupled many-body fermion systems. Signifi- 
cant progress in reproducing the empirical data has been made with the quantum 
Monte Carlo calculations for light nucleiu and with the variational chain summa- 
tions method for nucleon mattera. However, to ensure the proper binding energy 
of finite nuclei and the equilibrium properties of nuclear matter, these approaches 
make use of some plausible density-dependent or three-nucleon phenomenological 
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terms which are additionally included in the nuclear Hamiltonian. Then, with the 
present-day Monte Carlo calculations, an accuracy within 1% could be obtained for 
the binding energies of light p-shell nuclei up to the mass number A = 8, but, de- 
spite the rapid growth of computational power, for heavier nuclei such calculations 
are not feasible in nearest future. Moreover, unfortunately, this kind of "ab initio!' 
approaches (like the Brueckner-Bethe-Goldstone many-body theory, see e.g. Ref.B) 
do not yet provide an effective interaction, or nuclear energy-density functional, 
which could be used in nuclear structure calculations with sufficient accuracy to 
meet modern experiments and to give reliable predictions for nuclei far from the 
measured region. An important issue is how to reveal the density dependence of 
the effective interaction particularly in the particle-particle (pp) channel where the 
s-wave pairing correlations are known to play an essential role in finite nuclei. 

Among the existing effective models, the most successful are the phenomeno- 
logical mean-field microscopic approaches incorporating density-dependent forces 
of the Skyrme type with zero range or of the Gogny type with finite range, and also 
the relativistic mean field (RMF) model with classical meson fields. With some 10 
fitting parameters, these models can give the masses and radii of measured nuclei 
with the respective rms deviations ~ 2 MeV and « 0.02 fm from experimentQiJ. 
At the same time, they reveal a large spread in the extrapolation behavior to the 
unexplored regions of the nuclear chart and their predictions, already for nuclei 
not too far from stability, often deviate significantly from those of the macroscopic- 
microscopic (MM) mgdelsa or of the extended Thomas- Fermi model with Strutin- 
sky integral (ETFSI)EI. These latter models are able to reproduce the masses and 
charge radii of known nuclei with rms errors down to rs 0.6 MeV and w 0.02 fm, 
respectively. One may still notice that these errors are not mutually consistent: the 
relative rms deviations for radii are a few times larger than those for masses. Now, 
the description of the bulk nuclear properties within a microscopic model with the 
accuracy on the level of about 0.1% would be considered as quite successful. 

2 A new form of the LEDF 

Recently, a new ansatz for the LEDF construction has been attempted!! to diminish 
the gap between the predictions of the fully self-consistent microscopic mean-field 
models and the quasi-classical or MM models, hunting also for a more universal 
nuclear density functional which could be used not only throughout the nuclear 
chart but also for describing such objects as neutron stars, with crystal structure 
in their crust. The energy density of a nuclear system is represented as 

£ = £kin + £v + £s + ECoul + £sl + Eanom , (1) 

where e^in is the kinetic energy term which, since we are constructing a Kohn-Sham 
type functional, is taken with the free operator t = p 2 /2m, i.e. with the effective 
mass 77i* = 777,; all the other terms are discussed below. 
The volume term in ([l|) is chosen to be in the form 
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Here and in the following x± — (p n ± p p )/2po with p n (p) the neutron (proton) 
density, 2po is the equilibrium density of symmetric nuclear matter with e F = 
(97r/8) 2 / 3 ?i 2 /2mrg, the Fermi energy and r = (5/8-Kpo) 1 ^ 3 , the radius parameter. 
The fractional expressions of the type of Eq. (ph were introduced in Ref.El for the 
LEDF with application to finite systems with s-wave pairing correlations. Such 
expressions allow an extrapolation of the nuclear equation of state (EOS) to very 
high densities while preserving causal behavior. This might be of advantage since 
the available microscopic EOS often violate causality at p > 1 £m_ , 3 . Thus, in 
deriving the parameters of Eq. (||), we shall use the EOS of Refs.EjO only in the 
region of up to about six times the saturation density. The four parameters in the 
isoscalar volume energy .density tx a^_ are found by fitting the EOS of symmetric 
infinite nuclear matterErEJ for the UV14 plus TNI model. The result shown in Fig. [l] 
by the lower solid curve is obtained with the exponent a = 1/3, the compression 
modulus A'o = 220 MeV, the equilibrium density 2po = 0.16 fm~ 3 (tq = 1.143 fm) 
and the chemical potential p = —16.0 MeV (the energy per nucleon at saturation 
point). The corresponding dimensionless parameters are a+ — —9.559, h\ + — 
0.633, h\ + =j_D_L31. Keeping them fixed, a fit to the neutron matter EOS from the 
same paperslljO is performed to determine the three parameters of the isovector 
part oca^l in Eq. (@). The result shown in Fig. [H by the upper solid curve is obtained 
with al_ = 4.428, h\_ = 0.250, h v 2 _ = 0.130. This corresponds to the asymmetry 
energy coefficient (3q = 30.0 MeV. The surface part e B in Eq. ([!]) is meant to describe 
the finite-range and nonlocal in-medium effects which, phenomenologically, may 
presumably be incorporated within the LEDF framework in a localized form by 
introducing a dependence on density gradients. It is taken as follows: 
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Figure 1. The energy per nucleon versus density in nuclear and neutron matter. Open circles and 
crosses are the calculations by Friedman and Pandharipande 10 and Wiringa et oi. , respectively. 
The solid curves are obtained from Eq. fch with the parameters given in the text. 
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with h s , = ft-2+: an< ^ the two free parameters a s + and /i v . This peculiar surface 
term may be regarded as the Pade approximant for the (unknown) expansion in 
(Vp) 2 /(1 + h s + x+) where the forrn_factor 1/(1 + h s + x°,) imitates a transformation 
to Migdal's quasiparticles (cf. Ref.Ej). In fact, h s + is also a free parameter but here 
we prefer to keep it fixed by the above condition. 
The Coulomb part in Eq. (|lj) is approximated by 

£Coul = ^p ch (r) J g|g^ - ^ )1/3e2 ^ /3(l _ ^ oui<) ) (4) 

where the first term is the direct Coulomb contribution expressed through charge 
density p c h while the second term is the exchange part taken in the Slater approx- 
imation and combined with the Coulomb- nuclear correlation term ex /icoui- The 
latter is believed to account for the correlated motion of protons, in nuclei beyond 
the direct (Hartree) and exchange (Fock) Coulomb interactionEJ. The parameter 
/iCoui allows to practically kill one more enemy: the Nolen-Schiffer anomaly. 

The spin-orbit term e B i in Eq. (fil) is related to the two-body spin-orbit interac- 
tion oc (k + k'ti ■ f 2 ) [V 1 8{f\ — r 2 ) x (pi — p 2 )] • (<?i + 1?2 ) ■ It is known from the RMF 
theory that the isovector spin-orbit force is very small compared to the isoscalar 
one. Thus we set k! = and derive the isoscalar strength n = 0.19 from the average 
description of the splitting of the single-particle states in 208 Pb. 

The last term in Eq. (TO), the anomalous energy density, is represented as 

Eanom = £ V* (?)!* (x+ (f))lS (f) , (5) 

i— n, p 

where v(r) is the anomalous density and T^ = Cof^ is the^effective force in the 
particle-particle channel with the dimensionless form factortJ 

fHx+) = fL + h^xl + fir*(Vx+) 2 . (6) 

Here Co = 2ep/3po is the inverse density of states on the Fermi surface (Co = 
307.2 MeV-fm 3 ); q = 1. The strength parameters /| x = -2.8, hfi = 2.8 and 
/y = 2.2 are extractfid from a fit to the neutron separation energies and charge 
radii of lead isotopesliJ. 

The three parameters a s + , /i v and /icoui remain to be defined. This was done 
through a x 2 fit to the masses and radii of about 100 spherical nuclei from 38 Ca to 
220 Th with the result a s + = 0.600, h s v = 0.440 and h Co ui = 0.941, the rms deviations 
being 1.2 MeV and 0.01 fm for masses and radii, respectively (the zero-point energy 
correction — |41A -1 / 3 MeV has been included). rThe LEDF in the suggested form 
with the above parameters has been abbreviated^ as FaNDF . 

Typical results of the spherical LEDF calculations with FaNDF are shown in 
Fig. for even Pb and Sn isotopes, from the proton to the neutron drip line, in 
comparison with experimental data and other model predictions. The ETFSI model 
is chosen as a reference. It is seen that the predictions obtained with the Gogny force 
just outside the measured regions are in strong disagreement with other models; 
the RMF calculationsEj yield larger oscillations in the masses around experimental 
values and seem-.to predict different behavior in the neutron-rich domain compared 
to others (Rcf.Ej provides no results for the near-drip-line isotopes above A = 252 



ledfans: submitted to World Scientific on March 29, 2000 



14 
12 

> 10 
u 8 



E 2 

§ -4 
-6 



o exp. 

MM 


Pb 






Gogny 


/ 




/ 


-— RMF 


/ 




/ 


7 FaNDF 


/ /^-^ 




/ '. 


r / 


L^/^ 


r7 


■ 


*• «aa- , fai^^^^ e ^M 






■ 


f \M" 






■ 


/ 






- 


; / 

,<.. \ .... \ .... \ .... \ . 


, , , i . , , , i , , , , i . , 


. i , , , . i 


- 



170 180 190 200 210 220 230 240 250 260 270 



5 2 



; 



Sn 



exp. 
MM 



f \ 

y"~^, I ■ ■ Gogny 

,.- V ' / \ "— RMF 

1 FaNDF" 




100 110 120 130 140 150 160 170 
A 



Figure 2. Deviations of various theoretical masses from the ETFSf mass 7 for lead and tin isotopes. 
Dotted lines correspond to the MM model 6 , dashed to the Gogny force 5 and dash-dotted to the 
RMF model 15 . The results obtained with FaNDF are shown by solid lines. The experimentally 
known masses 16 are presented by open circles. 



for lead and above A — 156 for tin). Approaching the neutron drip line, the FaNDF 
masses fall in between the MM and ETFSI predictions. 

The role of the term oc /icoui entering Eq. (0) can be illustrated by the following 
example. The mass differences for the mirror pairs 17 F- 17 and 41 Sc- 41 Ca calcu- 
lated with FaNDF are 3.546 MeV and 7.174 MeV, respectively, whereas the corre- 
sponding experimental values are 3.543 MeV and 7.278 MeV. If one sets /icoui = 0, 
the calculated mass differences for these mirror pairs would be respectively 3.300 
MeV and 6.872 MeV leading to a 6-7% discrepancy known as the Nolen-Schiffcr 
anomaly. It follows that Coulomb-nuclear correlations play an important role in 
finite nuclei. Incorporating the corresponding term in the LEDF improves the de- 
scription of nuclear ground state properties and greatly reduces the severity of the 
Nolen-Schiffer anomaly. 
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3 Variational principle and coordinate-space technique 

The LEDF calculations discussed above were performed by using the diagonal ap- 
proximation in the pairing channel, with the particle continuum approximated by 
the quasistationary states, with each wave function being normalized in a finite 
volume which includes only the decaying tail beyond the nuclear surface. Such a 
HF+BCS technique eliminates to a great extent the presence of the "particle gas" 
and allows faster large-scale calculations. In certain situations, however, a more 
accurate treatment of the coupling with the particle continuum is needed (strongly 
surface-peaked pairing, nuclei near the neutron drip line, fine peculiarities in the 
behavior of differential observables, e.g. the staggering in charge radii, etc.). An 
appropriate tool is then the coordinate-space Gor'kov (or HFB) equations which 
naturally yield the localized ground state wave functions fr>r fi-nitn systems, with 
correct asymptotics for the normal and anomalous densitiestlfESJ. A more rigor- 
ous formulation of the LEDF approach based on the general variational principle 
and the coordinate-space technique is also of relevanceEj. In brief, we proceed as 
follows. 

Generally, the energy E of a nucleus with pairing is a functional of the general- 
ized density matrix R which contains both a normal component p and anomalous 
component v: 

E[R]=E kin [p}+E int [p,i>}, (7) 

where Skin [p] = Tr(ip) , and E- mt [p, v\ = ^int (normal) [p] + -^anomai [p, i>\ ■ The anoma- 
lous energy -Eanomai is chosen such that it vanishes in the limit v — ■» 0. In the weak 
pairing approximation |Af/cf| <C 1, which is usually the case for nuclear systems, 
one needs to retain only the first-order term ~ v 2 : 

£ an omal[p, i>] = \ (v ] H P W) = \ (*> + A) , (8) 

where J- pp is the antisymmetrized effective interaction in the pp channel and the 
parentheses imply integration over all variables. To calculate the ground state 
properties, one can now use the general variational principle with two constraints, 
(HFB|A r (^)|HFB) = N(p) = N and R 2 = R, leading to the variational functional 
of the form 

I[R] =E[R] -pN(p)-TrA(R-R 2 ), (9) 

where N is the particle number, u_ the chemical potential, and A the matrix of 
Lagrange parameters (see e.g. RefEil). In general, the gap equation is nonlocal 
and its in^inedium solution poses serious problems. We use the renormalization 
procedures by introducing an arbitrary cutoff e c in the energy space, but such 
that e c > ep; and by splitting the generalized density matrix into two parts, R — 
R c + S C R, where 8 C R is related to the integration over energies |e| > e c . The gap 
equation is renormalized to yield 

A = \H»c , (10) 
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where v c is the cutoff anomalous density matrix and .F| is the effective antisym- 
metrized pairing interaction in which the contribution coming from the energy re- 
gion |e| > e c is included by renormalization. With the Green's function formalism, 
for homogeneous infinite matter, it is shown that the variational function E — fiN 
does not change in first order in A 2 upon variation with respect to S C R. The total 
energy of the system and the chemical potential also remain the same if one imposes 
the above constrainlJV(^) = N for the cutoff functional. To a good approximation, 
as discussed in Ref.E3, this should be also valid for finite (heavy) nuclei. This result 
is expected since the major pairing effects are developed near the Fermi surface 
and the pairing energy is defined by a sum concentrated near the Fermi surface. In 
infinite matter the pairing energy per particle is E pa i T /N = — 3Ap/8ep. It follows 
that, with the cutoff functional, provided e c > ep, this leading pairing contribution 
is exactly accounted for. Thus the nuclear ground state properties can be described 
by applying the general variational principle to minimize the cutoff functional which 
has exactly the same form as (g) but with R replaced by R c . Recalling now the 
Hohenberg-Kohn existence theorem, one can choose the energy functional to be of 
a local form, the LEDF, dependent on the normal and anomalous local densities 
p(r) and v{f). Then, after separating the spin variables, the anomalous energy 
density acquires the form of Eq. (H), and (10) becomes simply the multiplicative 
gap equation 

A(r)=^(r;[p c ]K(r). (11) 

Having found the pairing and the HF potentials, tha-Gor'kov equations may be 
solved exactly with the coordinate-space techniqueEaEj. The Green's functions 
obtained this way are integrated in the complex energy plane to yield both the 
normal and anomalous densities which are used then to compute the energy of the 
system. It is worthwhile to notice that our approach does not imply a cutoff of 
the basis since the general variational principle is formulated with a "cutoff" LEDF 
from which the ground state characteristics of a superfluid system may be calculated 
by using the generalized Green's function expressed through the solutions of the 
Bogolyubov equations at the stationary point. To construct the densities that 
appear in this local functional, only those Bogolyubov solutions from the whole 
set are needed which correspond to the eigenenergies E a of the HFB Hamiltonian 
(which is a matrix of the first variational derivatives of the LEDF) up to the cutoff 
e c > e F . 

From the Gor'kov equations, after separating angular variables, for the gener- 
alized radial Green's function gji one gets the equation 

e-hji+n -A V I \ ( 5 (n-r 2 ) \ . . 

-A e + hjl -^)^(n,r2;e)=^ Q S{ri _ r2) ), (12) 

where hji is the single-quasiparticle HF Hamiltonian in the jl channel. The solu- 
tion of this matrix equation can be constructed by using the set of four linearly 
independent solutions which satisfy the homogeneous (Bogolyubov) system of equa- 
tions obtained from (O) by setting the right hand side to zero and which_abey the 
physical boundary conditions both for the bound and scattering statesfL2lEl 
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4 Nuclear isotope shifts 

The formulated approach has been applied to the combined analysis of the differen- 
tial observables such as neutron separation energies and odd-even effects in charge 
radii along isotopic chains. Reproducing the observed changes of geometrical char- 
acteristics of nuclei, first of all the staggering and kinks in charge radii, could shed 
light on the density dependence of the effective pairing interactiorc 3 ! This point 
may be illustrated by considering how the density changes when the pairing gap 
appears in nuclear matter. To leading order, one finds the following expression for 
the energy per nucleon near the saturation point (henceforth, x = x + ): 

A A 18 Xq 8 €f(x) 

where Kq is the compression modulus and (3(x)I 2 is the asymmetry energy with 
I = (p n — p p )/2pox = (N — Z)/A. Due to the pairing term oc A 2 , the position 
of the equilibrium point may be shifted to lower or higher densities depending on 
the behavior of A (a;) near x = 1. If A does not depend on x in the vicinity of 
x = 1 then the equilibrium density decreases due to presence of ep(x) ex x 2 ' 3 in 
the denominator of the pairing term (the system gains more binding energy). This 
means an expansion of the system. The effect is enhanced if A becomes larger 
during such an expansion, i.e. when the derivative dA{x)/dx at x = 1 is negative. 
At equilibrium the pressure P = -§-{E/A) = 0. For saturated symmetric nuclear 
matter without pairing the dimensionless density is, by definition, xq = 1. If / =/= 
and A^O, the new equilibrium density p = 2pox can be found from the equation 



too 

X - Xo = -r^~ 



A(x) f z dA(x) A(x)\ df3(x) l2 



4,€f(x) \ dx x ) dx 



(14) 



One can see that the density should be sensitive indeed to the derivatives of A, 
and a negative slope in A should cause a decrease of the density. The obtained 
relations can_be used to estimate the influence of pairing on the charge radii for 
heavy nucleiEx It was shown that pairing interaction with strong p-dependence at 
x rs 1 might significantly change the equilibrium density. The size of this effect is 
controlled by the parameter ht. For reasonable values of h£, as given by Eq. (n5f) 
below, the shift of the saturation point is relatively small: \Sp\/2po < 0.8 %. In 
finite nuclei, the surface term ex /y is equally important to produce a kink at magic 
neutron numbers and, especially, to explain the observed odd-even staggering in 
(f 2 )ch along isotopic chainalil. The density dependence of the pairing force leads to 
the direct coupling between the neutron anomalous density u n and the proton mean 
field. The suppression of v n in the odd neutron subsystem because of the blocking 
effect influences the proton potential through the volume, cx/i^, and surface, oc/|,, 
couplings, and this moves the behavior of the proton radii towards the desired 
regimeE-i The results for the lead chain obtained with the elaborated coordinate- 
space technique and with the LEDF parametrization DF3 from Ref.Ej are shown in 
Fig. pi The following parameter sets of the pairing force <M) are deduced (for the 
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Figure 3. Upper panel: neutron separation energies S n for lead isotopes. Lower panel: differences 
of mean squared charge radii S(r 2 ) c h with respect to 208 Pb as reference nucleus; 72% of the 
corresponding liquid drop values (using rrj = 1.1 fm) are subtracted to enhance the visibility 
of the small differences. Curves (a)— (f) correspond to the parameter sets (a)-(f) of the pairing 
force ([19), respectively. Experimental data for S n are from Ref. 16 , for <5(r 2 ) c h from Refs. 26-28 . 
For comparison, the HFB calculations with Skyrme forces SkP and SLy4, and with the Gogny 
D1S force but for even isotopes and 211 Pb only, are also shown. 
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(15) 
As seen in the upper panel of Fig. S, the neutron separation energies ^n are repro- 
duced equally well for all parameter sets (a)-(f ) of Eq. (H) . The HFB calculations 
with two state-of-the-art Skyrme functionals SkP (RefTS) and SLy4 (Ref.cJ) de- 
scribe the S n values with more or less the same quality as our curves (a)-(f), some 
deviations are observed above 207 Pb. In the lower panel the calculated isotope shifts 
of mean squared charge radii, <5(r 2 ) c h, with respect to 208 Pb as a reference nucleus 
are presented. Our curves except (a) do not differ from each other significantly and 
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all of them reproduce qualitatively the kink and the size of odd-even staggering. 
The curve (a) corresponds to a simple contact pairing force which does not depend 
on density; in this case the behavior of <5(r 2 ) c h is rather smooth and neither kink 
nor staggering are reproduced. One can see also that the Skyrme functionals SkP 
and SLy4 give too small kink and too small staggering. The Gogny D1S force could 
not reproduce the kink either. This fact clearly demonstrates the importance of 
including odd-mass nuclei in the fitting procedure of the force constants. 

5 Ground-state properties of nuclear matter 

Now, the empirical information gained from finite laboratory nuclei can be used 
to study the ground state properties and behavior of A as a function of density 
p = 2pox (or the Fermi momentum k F = (3ir 2 p/2) 1 ' 3 = fcoF^ ) m infinite nu- 
clear matter. Since the density-gradient term oc /y vanishes in this case, only two 



parameters are relevant: /| x and h>. The gap equation (11) reduces to 



A(x) = -[ T^Jtjx) - A(X) , (16) 

A< fec (2tt) 3 J 2^(e k -e F (x)) 2 +A*(x)' 



where k c = y/2m(e F + cAJh and e k — U k 2 /2m. Its solution in the weak pairing 

regime can be written asE2l: 

r 7r i 

A(fc F ) = cepexp — — cotS(k F ) , (17) 

where c = 8e -2 s=s 1.083 and where we have introduced the Fermi level phase shift 
S(k F ) defined by 

, ,c/, x 4fc 0F / 1 , k c (k F )\ k F fk c (k F )-k F \ 

with k c (k F ) — \/fcg c + kp] fcoc = \ / 2me c /h. Eq. (118) corresponds to an exact solu- 
tion of the nn scattering problem at the relative momentum k — k F with the states 
truncated by a_momentum cutoff fc c = k c (k F ) for contact potential Cof^(k F )5(r) 
(see, e.g., Ref.0). 

Shown in Fig. [| are the results for A in nuclear matter with parameter sets of 
Eq. ([l|). Interestingly, for the sets (b)-(e), which reproduce satisfactorily both the 
S n and <5(r 2 ) c h values, there exists a "pivoting" point at k F rs 1.15 fm" 1 (at «0.65 
of the equilibrium density) with the same value of A p i v s=s 3.3 MeV. The approxi- 
mation ( |l7| ) (dashed curves in Fig. [I|) works well in the entire range of k F for the set 
(a), but for the other sets this is true only at k F greater than ssl.2 fm _1 and also, 
for the sets (b), (c) and (d), at k F less than 0.42, 0.14 and 0.042 fm -1 , respectively 
(in these regions, A/ep < 0.1). It should be stressed that e F entering the integrand 
of Eq. ( |16| ) can be expressed directly through density p by e F = h 2 k F (p)/2m with 
k F {p) = (37r 2 p/2) 1 / 3 only if the pairing is weak and the dependence of the Fermi 
energy e F (and the chemical potential p) on A can be disregarded. Otherwise one 
should introduce the particle number condition 

2 f d ^ ( \ (~\ 1 1 1 e k -e F (x) \ 

Tn k (x), n k (x) = - 1 = = == = , (19) 



PO J k < kc (2^) 3 W ' 2 V ^{e k -e F {x)Y+^{x) 
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Figure 4. Pairing gap in nuclear matter as a function of the Fermi momentum. Full (dashed) curves 
(a)-(f)_are obt ain ed by solving Eqs. ( Jlq ) and i 191) (by using the weak pairing approximation of 
Eqs. jl7[ ) and (D_q)) with contact pairing force (pjand energy cutoff e c = 40 MeV, and correspond, 
respectively, to the parameter sets (a)— (f) of Eq. (D-9). The vertical arrow marks the Fermi 
momentum at saturation point for the functional DF3 without pairing, fcoF = 1-328 fm — *. The 
solid (open) circles are the solutions of the nonlocal gap equation with the CD-Bonn potential 
(with the finite-range Gogny Dl force 33 ). The dotted line is obtained from Eq. JIT] ) with the free 
NN scattering phase shifts (see text). 




1.5*10 



P (fm J ) 



Figure 5. Pairing gap A as a function of p at very low densities. Curves (c)— (e) correspond to the 
parameter sets (c)-(e) of Eq. (D_q), respectively. 



and to solve the system of the two equations (JLq) and (h9J) with respect to A and ep- 
The results shown in Fig. H by the solid curves correspond to such a solution. All 
parameter sets rtlq) except (a) reproduce the neutron separation energies and the 
isotope shifts of charge radii (r 2 ) cn of lead isotopes fairly well (see Fig. ||). Shown 
also in Fig. are the values of the 1 Sq pairing gap injiiiclcar matter obtained for 
the CD-Bonn potentiate and for the Gogny Dl forceH. The agreement between 
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the two latter calculations is relatively good, while both deviate noticeably from 
our predictions. The curve for density-independent force, set (a), stands by itself 
with a positive derivative dA(xl/dx everywhere; no acceptable description of (r 2 ) c h 
could be obtained in this caseo. 

At very low densities Jlq) reduces to 



&f cot 5(k-p) 



1 



""nn^F 



2fcp 

7T 



k F 2/i« 

2kVc ~ (fL) 2 



where a nn is the singlet scattering length, 



\/2me c 



2fc F V ^ & of /, 



4fc F V/< 



f 4 

/ox 



_fcp_\ 

koFj 



1 

./ cr 



3<j-l 



(20) 



(21) 



and r nn is the effective range, r nn = 4/7rfc uc . Here we have introduced the criti- 
cal constant /| r = — 2fcoF/^Oo the vacuum strength /| x at which the two-nucleon 
problem has a bound state solution at zero energy (in our case, /| r = —1.912). 

The first two terms in ( p0|) would describe the s-wave phase shift at low energies 
through an expansion of k cot S in powers of the relative momentum k = kp if the 
interaction were density-independent - i.e., if the coupling strength and momentum 
cutoff were fixed by /* = /| x and k c = fco c , respectively. It follows that, with 
a density-dependent effective force, such an expansion contains additional terms 
which are, for the parametrization used here, of the same order as the effective 
range term. This simply demonstrates that, for reproducing the pairing gap, the 
effective interaction even at very low densities need not necessarily_cpincide with 
the bare NN interaction as was discussed by Migdal many years agco. 

At very low densities, at fcp — ► 0, to leading order from (O) we obtain 



ceF exp 



a nn < . 



(22) 



This expression agrees with the results of Ref£3 based on a general analysis of the 
gap equation at low densities^] when kp \ a nn | < 1. 

But we should stress that (|22j ) is valid only in the weak-coupling regime cor- 
responding to negative a nn . In the opposite case, the gap in the dilute limit has 
to be found in a different way. At /| x > /| r , from (J2^) and (|l^) it follows that at 
low densities the pairing gap is exponentially small and eventually A(fcp — ► 0) = 0. 
Such a weak pairing regime with Cooper pairs forming in a spin singlet I = state 
exists up to the critical point at which the attraction becomes strong enough to 
change the sign of the scattering length. Then the strong pairing regime sets in, and 
A should be determined directly from the combined solution of Eqs. ( |l6| ) and (|l9|). 
In the dilute systems, €f plays the role of the chemical potential \i. At the critical 
point, /i becomes negative, and a bound state of a single pair of nucleons with the 



a Our effective pairing interaction with the choice q = 1/3 would lead in the dilute limit to the 
expression for A of the form of Eq. (E3) but with a different prefactor c depending on the value 
of h^. If, furthermore, we define it by h^ = (1 + 21n2)(/,J x )fV|6 we get in the leading order 
A(fcp) = (2/e) 7 ' 3 ep exp(7r/2fcpa), i.e. the result obtained in RefE3 for a non-ideal Fermi gas with 
taking into account the terms up to the second order in kpa. 
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binding energy e& = 2/x (= —h/ma^ n ) becomes possibleEj. In this regime, A can 
be found from (191). In the leading order we get 



A = -(^l , a nn >0. (23) 

m \ a nn J 

It follows that, in the dilute case, the energy needed to break a condensed pair goes 
smoothly from 2A to 2p = tb as a function of the coupling strength as the regime 



changes from weak to strong pairing. But as seen from (22) and (G3J), the behavior 
of A is such that the derivative dA/dp at p — > as a function of /| x exhibits a 
discontinuity from to oo. This is illustrated in Fig. where we have plotted 
A(p) for the sets (c)-(e) of Eq. (|l5|), which embrace both regimes. Notice also that 
both analytical expressions (E2I) and ( p3| ) give a pure imaginary gap at the critical 
point, where the scattering length changes sign. When the Fermi momentum k F 
approaches from below the peper critical point, where the pairing gap closes, A 
becomes exponentially smallti In weak coupling, A is also exponentially small at 
low densities. It is noteworthy that, as follows from (J17|), in both these cases the gap 
closes exactly at the points where the phase shift passes zero. As an illustration, we 
show in Fig. [| by the dotted line the values of A(fcp) obtained from ( |17|) using the 
"experimental" nn phase shiftsg. It is seen that A obtained ±his way closely follows 
the solution of the gap equation with the CD-Bonn potential^ 2 ] at low densities. The 
nn phase shift passes zero at the relative momentum k = 1.71 fm _1 , and the gap 
should vanish at the corresponding Fermi momentum. Unfortunately, the solutions 
for A are given in RefE3 only in the region up to k F = 1.4 fm _1 . For nuclear 
matter, with our LEDF, the energy per particle is 

E , n 2 /• dk h 2 k 2 . . 1 v „, , 3A 2 (ir) 

A pax J k < kc (27r) J 2m 3 2fi{x)xe QF 

where fl(x) = (1 -=Jf[ + x) / (1 + h% + x) . Numerically, a v + = -6.422, h\ + = 0.163 and 
h\ + = 0.724 (Ref.E3). In the dilute limit the "particle-hole" term ex /+ vanishes 
linearly in density. The chemical potential is 

1 3f*'M A 2 M 

»{x) = e F (x) + -e 0F al[fl'(x)x 2 + 2fl(x)x] + T ) l ] , (25) 

•J Z J (x) £of 

where the prime denotes a derivative with respect to x, e F (x) and A(x) are deter- 



mined from (16) and (|19|). The two last terms in (|25|), even in strong pairing regime, 
vanish in the dilute limit at least as x q if < q < 1 or linearly in x if q > 1 . We see 
again that, in strong coupling, in the leading order fi — e F = e^/2 < 0. The calcu- 
lated energy per nucleon as a function of the isoscalar density p is shown in.the upper 
panel in Fig. o together with the results of the nuclear matter calculationsc2lO. It is 
seen that DF3 gives qualitatively reasonable description of the nuclear matter EOS 
and that pairing could contribute noticeably to the binding energy especially at 
lower densities. In the lower panel in Fig. B we have plotted the pairing energy per 
nucleon, (E/A) pa ,i T , obtained by subtracting from ( pi} ) the corresponding value of 
E/A at A = 0. The pairing contribution increases, as expected, when /| x becomes 



We thank Rupert Machleidt for providing us with these nn phase shifts. 
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gradually more attractive, with a shift to lower densities. For the sets (e) and (f) 
the attraction is strong, /| x < /| r . In these cases a nonvanishing binding energy 
in the dilute limit is solely due to Bose-Einstein condensation of the bound pairs, 
the spin-zero bosons, when all the three quantities, /i, E/A and {E / A) pa \ Vl reach 
the same value eb/2 = fr 2 /2ma„ n (ej, = —0.0646 and —1.616 MeV for the set (e) 
and (f), respectively). This is illustrated in Fig. ^ where we have plotted E/A and 
(E/A) pa i T as functions of p at very low densities. 

We have considered the properties of nuclear matter with s-wave pairing within 
the LEDF framework, including extrapolation to the dilute limit, with a few possible 
parameter sets of the pairing force deduced from finite nuclei. At low densities, in 
the case of symmetric N = Z matter, the 3 S\ — 3 Di pairing could be more important 
since the n-p force is more attractive than in the p-p or n-n pairing channels. 
Thus, our approach, with the 1 Sq pairing only, would be more appropriate for an 
asymmetric N ^ Z case and for pure neutron systems. From this point of view the 
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Figure 6. Energy per nucleon E/A (top) and pairing contribution to E/A (bottom) in nuclear 
matter. Curve s (a .)— (e) are calculated using Eq. (E4) and correspond to the strength parameters 
(a)— (e) of Eq. (|15[), respectively. Open circles and stars are the calculations of Ref. 10 and Ref. 11 , 
respectively, for the UV14 plus TNI model. The vertical arrows mark the saturation density for 
the functional DF3 without pairing, 2po = 0.1582 fm — 3 . 
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Figure 7. Energy per nucleon E/A (top) and pairing contribution to E/A (bottom) for nuclear 
matter at low densities. The notations are the same as in Fig. 0. 



best choice for the LEDF calculations seems to be the parameter set (d) of Eq. ( p_5[ ) 
since it gives the singlet scattering length a nn rs —17.2 fm which corresponds to a 
virtual state at rs 140 keV known experimentally. As seen in Fig. [|, for this choice 
the behaviour of A at low densities agrees well with the calculations based on 
realistic NN forces. At higher densities, however, our predictions for A with the set 
(d) go much higher reaching a maximum of « 4.84 MeV at kp sw 0.92 fm -1 , while the 
calculations of Ref.ES give a maximum of about 3 MeV at kp « 0.82 fm -1 . With 
a bare NN interaction, if one assumes charge independence and that m n = m p , 
the pairing gap would be, at a given fcp, exactly the sam£-.lpath in symmetric 
nuclear matter and in neutron matter. As shown in Refs£§E!l, if one includes 
medium effects in the effective pairing interaction, A in neutron matter would be 
reduced substantially, to values of the order of 1 MeV at the most. Whether such 
a mechanism works in the same direction for symmetric nuclear matter is still an 
open question. The force (@) was chosen to be dependent on the isoscalar density 
only, since we have analyzed the existing data on S n and <5(r 2 ) c h for finite nuclei 
with a relatively small asymmetry, (N — Z)/A < 0.25. An extrapolation to neutron 
matter with such a simple force would give a larger pairing gap than for nuclear 
matter. This suggests that some additional dependence on the isovector density 
p n — p p might be present in the effective pairing interaction. 
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6 Staggering and kinks in nuclear multipole excitations 

The staggering phenomenon observed in the behavior of charge radii plotted as a 
function of neutron number is a prominent odd-even effect which has been system- 
atically measured and widely discussed. Similar effects are expected to exist for 
other quantities such as neutron and matter radii, centroid energies of multipole 
excitations (position of the giant resonances), etc. Shown in the upper panel of 
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Figure 8. Neutron, proton and matter radii (top), mean energy of isovector dipole (middle) and 
Fermi charge-exchange (bottom) transitions for tin isotopes. The solid dots for aJdip correspond 
to the experimental data for the maximum of the dipole photoabsorption Lorentz curves 44 , while 
those for Hip to the experimental positions of the isobaric analog states with respect to the daughter 
nuclei 45 . 
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Fig. |8| are the rms proton, neutron and matter radii for the tin isotopic chain cal- 
culated with the functional DF3 and with the set (d) of the pairing force. In the 
middle panel we show the mean energies <I>di P of dipole isovector excitations. These 
are calculated as the square root of the ratio m^/mi with mi and 7713 the linear 
and cubic energy-weighted sum rule, i.e. the first and the third moment of the 
corresponding RPA strength distribution, respectiveljo: 



Wdip 



3m NZ 



/*"*-' %^>-')f? 



1/2 

(26) 



Here J- np is the effective neutron-proton interaction obtained from the energy- 
density functional as the second variational derivative 5 2 E int /Sp n Sp p . In the lower 
panel of Fig. H we show the mean energies uj-p of the charge- exchange + excitations, 
i.e. the Fermi transitions, in tin nuclei with N > Z. These energies are calculated 
within the sum rule approach as the ratio (m{ + ra^)/{m^ — tUq) with mj and 
m7/ the first moment of the strength distribution of the Fermi transitions in the 
(3 + and /3~ channel, respectively (energy-weighted sum rules) and with m[j~ and 
tUq the corresponding non-energy weighted sum rules. The expression for uip reads 
(see, e.g. RefMd): 

^f = -j^ j dfU Coul (r) ( Pn (r) - Pp (r )) , (27) 

where C/coui is the Coulomb mean field potential. It is seen in Fig. g that the 
neutron, proton and matter radii as functions of the mass number A reveal a kink 
at magic 132 Sn, and at larger A the difference between rms neutron and proton 
radii starts to increase more rapidly. The staggering in radii is observed mostly 
in the region between the two magic nuclei, from 100 Sn to 132 Sn, and this effect 
is practically washed out beyond A = 140. The mean energy of dipole transitions 
occurs to be in anticorrelation with such a behavior in radii, and this seem to be 
in qualitative agreement with experimental data (one should mention that Eq. (Eq) 
overestimates the position of the giant dipole resonance by rs 0.5 MeV because 
of the cubic energy- weighted sum rules used in the derivation of u>dip)- One also 
observes a distinct kink in the behaviour of a>di P at A = 132. Beyond this magic 
number the mean dipole energy decreases rather fast and then nearly saturates 
when approaching A = 172. Ihis might be connected with an enhancement of the 
low-energy dipole-.transitionsc 2 ] and also with possible appearance of the so-called 
soft dipole modeCJ in nuclei near the neutron drip line. The anticorrelations in the 
staggering behavior of tZ>di P and the rms radii (r 2 ) n , (r 2 ) p can easily be under- 
stood by considering the influence of pairing on the gradients of neutron and proton 



densities entering Eq. (26). The odd-even effect in <Ddi P is constructive and more 
pronounced in the A < 132 region where both gradients strongly overlap. Beyond 
A = 132 the larger differences between neutron and proton rms radii imply a lower 
overlap between neutron and proton density gradients at the nuclear surface, hence 
a smaller mean dipole energy. The situation with mean energy of the Fermi charge- 
exchange transitions is different. From Eq. (p7J) one expects that the correlation 
between tip and the neutron and proton rms radii should be destructive. As seen 
in the lower panel in Fig. pL the staggering in the evolution of (Dp with A is very 
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weak indeed and almost invisible. The kink at the magic mass number A — 132 is 
also much less pronounced compared to the dipole case. Remarkable enough, the 
theoretical self-consistent sum-rule predictions for dip in tin isotopes are in excel- 
lent agreement with the available experimental data on the position of the isobaric 
analog states. 
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